Intrinsic alignment-lensing interference as a contaminant of cosmic shear 
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Cosmic shear surveys have great promise as tools for precision cosmology, but can be subject 
to systematic errors including intrinsic ellipticity correlations of the source galaxies. The intrinsic 
alignments are believed to be small for deep surveys, but this is based on intrinsic and lensing 
distortions being uncorrelated. Here we show that the gravitational lensing shear and intrinsic 
shear need not be independent: correlations between the tidal field and the intrinsic shear cause the 
intrinsic shear of nearby galaxies to be correlated with the gravitational shear acting on more distant 
galaxies. We estimate the magnitude of this effect for two simple intrinsic alignment models: one in 
which the galaxy ellipticity is linearly related to the tidal field, and one in which it is quadratic in 
the tidal field as suggested by tidal torque theory. The first model predicts a gravitational-intrinsic 
(G/) correlation that can be much greater than the intrinsic-intrinsic (//) correlation for broad 
redshift distributions, and that remains when galaxies pairs at similar redshifts are rejected. The 
second model, in its simplest form, predicts no gravitational-intrinsic correlation. In the first model 
and assuming a normalization consistent with recently claimed detections of intrinsic correlations 
we find that the GI correlation term can exceed the usual // term by > 1 order of magnitude and 
the intrinsic correlation induced B-mode by 2 orders of magnitude. These interference effects can 
suppress the lensing power spectrum for a single broad redshift bin by of order ~ 10% &t Zs = 1 
and ~ 30% at Zs = 0.5. We conclude that, depending on the intrinsic alignment model, the GI 
correlation may be the dominant contaminant of the lensing signal and can even affect cross-spectra 
between widely separated bins. We describe two ways to constrain this effect, one based on density- 
shear correlations and one based on scaling of the cross-correlation tomography signal with redshift. 

PACS numbers: 98.80.Es,98.62.Sb,98.62.Gq 



I. INTRODUCTION 

Weak gravitational lensing has attracted considerable 
interest recently as a means of measuring the density 
perturbations in the universe at low redshifts because 
it is directly sensitive to the matter distribution. It 
thus avoids the astrophysical complications involved in 
other means of measuring the matter distribution such 
as galaxy clustering. Recent first detections of the cos- 
mic shear autopower spectrum [ESIljSlljIaBjH have 
stimulated proposals for more ambitious pro ject s in the 
future, such as CFHTLS [H, Pan-STARRS 0, SNAP 
|6ll |. and LSST [62|, which aim for percent-level preci- 
sion. Comparison of the shear power spectra at differ- 
ent redshifts and/or to the primordial fluctuation ampli- 
tude derived from cosmic microwave background (CMB) 
anisotropy measurements, can constrain the cosmologi- 
cal growth factor and hence parameters such as the neu- 
trino mass, amplitude of fluctuations as, and dark energy 
equation of state H IH IH El El El . 

While the underlying physics of weak lensing is 
"clean," the shear measurements are subject to possi- 
ble systematic errors including incomplete correction for 
seeing and optical distortions, selection effects, and noise- 
rectification biases EE EE E3, EM El I^S I and their cos- 
mological interpretation relies on accurate knowledge of 
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the redshift distribution of the source galaxies. Another 
possible systematic error is intrinsic (i.e. not Icnsing- 
induced) correlations among the ellipticities of neighbor- 
ing source galaxies IHHllililiEililiiliill^, 
which could arise if the galaxy ellipticities are affected 
by large-scale tidal fields. This systematic is particu- 
larly worrisome because it lies outside the control of the 
observer, and is dependent upon the poorly understood 
physics of galaxy formation. One frequently proposed 
method to reduce this contamination is to assign photo- 
metric redshifts to the source galaxies, and then to down- 
weight or ignore pairs of galaxies at similar redshifts when 
computing the shear correlation function or power spec- 
trum ,31^ ^, ,3^ ^] ; see Ref . ;35|| for an implementation. 
The idea is that the pairs of galaxies widely separated 
in redshift should have independent intrinsic ellipticities 
but should have correlated gravitational shears induced 
by structures between the observer and the more nearby 
source galaxy. 

While cross-correlation of different redshift bins is ex- 
pected to remove spurious power due to the intrinsic 
alignment autocorrelation, there is another, more subtle 
effect by which even the cross-correlations can be con- 
taminated from intrinsic alignments. If the ellipticities 
of galaxies are correlated with the tidal quadrupole field 
in which the galaxies form, then the intrinsic ellipticity of 
a nearby source galaxy will be correlated with the lens- 
ing shear acting on a more distant source galaxy. This 
leads to a nonzero cross-correlation between the intrin- 
sic ellipticity and the gravitational lensing shear. Cross- 
correlations between the shear measurements at widely 



different redshifts are actually more contaminated by this 
effect than shear autocorrelations computed using only 
galaxies in a narrow redshift slice, because the radial sep- 
aration is necessary in order for the tidal field around 
the nearby galaxy to lens the more distant galaxy. This 
intrinsic-lensing correlation will also affect attempts to 
cross-correlate cosmic shear surveys of galaxies with lens- 
ing of the CMB 13, 36, 3^1 because the CMB is also 
lensed by the tidal field surrounding the source galaxies. 
This paper is organized as follows. In Sec.^ we for- 
mally express the E- and B-mode shear power spectra 
in terms of the background cosmology, power spectra 
of matter density and intrinsic shear, and matter-shear 
cross-spectrum. In Sec. Illll we consider two crude mod- 
els of intrinsic alignments and calculate their predicted 
contribution to the shear power spectrum. We discuss 
methods to assess and/or remove the contamination in 
Sec. IIVI and we conclude in Sec. |3 




FIG. 1: The effect of the density-intrinsic shear correlation 
on the shear power spectrum. Density fluctuations in the 
nearby plane (gray masses) induce a tidal field (arrows). A 
source galaxy in a more distant plane (dashed ellipse) is grav- 
itationally sheared tangentially to these masses. If the intrin- 
sic shears of galaxies in the nearby plane (solid ellipse) are 
aligned with the stretching axis of the tidal field, then this 
results in an anti-correlation between the shears of galaxies 
at different redshifts, i.e. C^ ' < 0. (The opposite case, 
Cf ' > 0, results if galaxies are preferentially aligned with 
the compressing axis of the tidal field.) 



II. SHEAR POWER SPECTRA 

Before developing the formalism let us describe a sim- 
ple example of the effect, shown in Fig. ^ The tidal 
field may lead to a stretching of the galaxy shape in the 
direction of the tidal field. Gravitational shearing of a 
background source leads to stretching of the galaxy a in 
perpendicular direction. As a result, the lensing effect 
will be partially cancelled by the intrinsic alignemnt ef- 
fect and the two effects are coherent, as they depend on 
the same underlying density field. 

Weak gravitational lensing by large-scale structure is 
detectable through its shearing of distant "source" galax- 
ies. To lowest order, the shear of a galaxy i can be broken 
down into a gravitational and an intrinsic shear contri- 
bution: 7i — 'ff + 7/. The gravitational shear is well 
known and is equal to: 



(7f+,7fx)=5- 



Wix,X^M-^l2^,^y)Sixn^)dx, 



where the d derivatives are takes with respect to angular 
position (i.e. have units of radians"^), d^^ is the as- 
sociated inverse Laplacian, hi is the angular position of 
galaxy i, S{x'0-) is the fractional density perturbation at 
distance x in direction n, and the lensing window func- 
tion is: 

3 
Wix,X^) = -n„^H^{l+z)sm% x{cotK X~cotK Xi) (2) 

for X < Xi 3'iid otherwise. Here sitik and cotK are the 
modified trigonometric functions, i.e. 

( K-'/^&m(K'/^x) K>0 

siiiK X^ { X K = , (3) 

[ |if|-i/2sinh(|if|i/2^) K <0 

cot if X = ^ Insiui^- X, and K is the spatial curvature of 
the universe. 



Now we imagine that a set of source galaxies in red- 
shift slice a with comoving distance distribution faix) 
are observed. The gravitational contribution to a shear 
Fourier mode is: 

/•CXD 

7p(a) = (cos2(/.i,sin20i) / W„{x)Si{x)dx, (4) 

Jo 

where the integrated window function is: 

/•OO 

M^o(x)= / fa{x')W{x,x')dx'. (5) 

Jo 



by 



The intrinsic alignment contribution at a point is given 



-f\h,a)^ I faixW [xi^)dx. 



(6) 



where the density- weighted intrinsic shear 7^ = {\-\-5g)^^ 
is computed from the fractional overdensity of galaxies 



(1) 5g and average intrinsic shear of galaxies 7^. A density 



weighting is technically necessary in Eq. (0) , which makes 
Wa{x) slightly dependent on angular position. On sub- 
arcminute scales where the fluctuations in 5g are large 
this results in production of i3-modes in the lensing shear 
[38j : on larger scales the effect is unimportant. In con- 
trast, intrinsically aligned pairs of galaxies tend to be 
close to each other where i5g > 1, and hence the factor 
oi \ + 5g in Eq. (jHl cannot be safely neglected except on 
very large scales. 

The -E-mode shear cross-spectrum between two red- 
shift slices can be broken down into gravitational lens- 
ing (GG), intrinsic alignment (//), and interference (GI) 
terms: 

Cf ^(a/3) = Cf ^'<=<=(a/3) + Cf ^'"(a/S) + Cf ^'^^(a/?). 

(7) 
(The _B-mode shear cross-spectrum is similar, but con- 
tains only an // term since there is no gravitational 



contribution to the J5-mode shear.) The gravitational 
lensing contribution to the shear power spectrum can be 
obtained on small scales {£ ^ 1) by Limber integration 
[39j, which computes the power spectrum of a radially 
projected quantity. In this case we integrate Eq. Q) to 
get 



Cf^'^^{al3) 






sm^X 



P&{^= 771 — -r:^x)dx- 

(8) 

The intrinsic alignment contribution is most easily com- 
puted by Limber integration of Eq. © : 






smJfX sm/fX 



where the projected power spectrum of the intrinsic 
alignments is: 

P^'^'ik-.X)^ n 1 1 {{^imiixWoikx^) + Mkx^)] 
+ {limiix))[Jo{kx^) - Mkx^)]} 



xxj_dxj_dxii, 



(10) 



where the () are two-point correlation functions of 7 with 
two points separated by comoving distance X|| in the ra- 
dial direction and x± in the perpendicular direction. The 
+ and X components of the shear are measured along the 
direction of separation x_l and at 45 degrees to this di- 
rection, respectively. Note that in general the two-point 
correlation functions can depend on x±_ and x\\ indepen- 
dently (i.e. they are not just functions of the separation 

\/x^j_ + X?. ) since the shear depends on viewing geometry. 

To get the i3-mode power, switch the -I- and x labels in 
Eq. Iini). 

Given models for intrinsic alignment and galaxy clus- 
tering, Eqs. (O and H10() are sufficient to compute the 
pure intrinsic alignment contribution to the shear cross- 
spectrum. However, note that as a second-order statistic, 
Eq. {Tj) also contains interference terms due to correla- 
tion of the gravitational shear with intrinsic alignments. 
These have normally been ignored because the gravita- 
tional shear acting on a particular galaxy is determined 
by the integrated tidal quadrupole along the line of sight, 
whereas the intrinsic shear is expected to be determined 
by tidal fields in the vicinity of the galaxy. But in cos- 
mic shear studies, we correlate the shears of two galaxies 
whose redshifts may be different. In this case, it is plau- 
sible that the more nearby galaxy is intrinsically aligned 
by a quadrupolar tidal field that also lenses the more dis- 
tant galaxy. Therefore we must consider the interference 
terms. These are given by 



Cf^'^\a[3) 



W^a(x)/Mx)+W^/3(x)/a(x) 

sin^ X 

xP5,^i(— ^)dx, (11) 

smx X 



where the 3-dimensional cross-spectrum Ps^=yi is related 
to the density-shear correlation function via 



-P<5,7j(fc;x) = -27r / / {5{0)^j^{x))J2{kx_[_)x_Ldx±dx\\ 



(12) 



III. MODELS 



We now compute the predicted level of contamination 
of weak lensing surveys from the "interference" terms, 
Eq. (|11() . We consider three cases: first the case of galax- 
ies whose mean ellipticities are linear functions of the 
underlying tidal field; and secondly, the case of galax- 
ies whose mean ellipticities are quadratic functions of 
the tidal field. In all of our numerical results below, we 
have taken the bias to be 6g = 1, used the best-fit scale- 
invariant (n^ = 1) flat ACDM model of Ref. [i^] {as = 
0.966, ^b = 0.0475, O^ = 0.293, Hq = 70.8 km/s/Mpc), 
and used the adiabatic CDM-|-baryon transfer function 
of Ref. 41] . The gravitational-gravitational shear power 
spectra include the nonlinear correction of Ref. [43 | . 

The first (linear) model has usually been used for el- 
liptical galaxies (since gravitational collapse in a gen- 
eral tidal field is expected to result in a triaxial halo 
aligned with the principal axes of the tidal field) and 
the second (quadratic) model for spirals (since the orien- 
tation of these is believed to be determined by angular 
momentum acquired during gravitational collapse, which 
requires one tidal quadrupole to supply the torque and 
another tidal quadrupole to be torqued). The models 
presented here are "toy" models in the sense that their 
theoretical motivation does not apply to highly nonlin- 
ear scales, and even on larger scales there is little ob- 
servational basis for models of intrinsic galaxy ellipticity 
alignments. An alternative to analytic models is to esti- 
mate intrinsic alignments by calculating the ellipticities 
or angular momenta of dark matter haloes in A^-body 
simulations [2l|, |23, 123, |23| , however one should keep in 
mind that there can be misalignment between the dark 
matter halo and the galaxy it contains [43, and so in- 
trinsic alignment results based on A'^-body simulations 
are not definitive. In principle it would be possible to 
derive predictions for intrinsic correlations from a halo 
model, however this would require adding a prescription 
for the galaxy alignment to the current halo models. Such 
a prescription would likely have to come from an analytic 
or A^-body model. Undoubtedly there is much room for 
improvement in both the theory - which must ultimately 
address the complications of galaxy formation, mergers, 
etc. - and the observations, which will be necessary to 
measure or constrain the intrinsic alignment signal at lev- 
els suitable for precision cosmology. 



A. Linear alignment model 

A simple model for the ellipticities of elliptical galaxies 
was proposed by [2g ■ The intrinsic shear of the galaxy 
is assumed to follow the linear relation 



/ = 



^' (V2-V5,2V,V,)5[*P], 



A-kG 



(13) 



where ^ p is the Newtonian potential at the time of 
galaxy formation (assumed to be during matter domi- 
nation), 5 is a smoothing filter that cuts off fluctuations 
on galactic scales, and V is a comoving derivative (as 
opposed to d which is a 2-dimensional derivative on the 
unit sphere with units of radians^^). We have taken <S to 
be a simple cutoff in Fourier space at kcutoS = 1/i/Mpc. 
Here Ci is a normalization constant (note that Ci > 
if the galaxy is aligned along the "stretching" axis of the 
tidal field). The original motivation for Eq. (|13|l was 
the assumption that galaxies are homologous with their 
haloes, and that the halo ellipticity is perturbed by the 
local tidal field produced by large scale structure |2^. 
On sufficiently large scales, this relation can also be moti- 
vated by arguments analogous to linear biasing theory for 
galaxies ^J: the large-scale correlations in the intrinsic 
shear field must be determined by the large-scale poten- 
tial fluctuations; if these large-scale potential fluctuations 
are sufficiently small, then the intrinsic shear field should 
be a linear and local function of the potential ^ p (which 
is oc ^P in the linear regime); the only linear, local func- 
tions of \l/p with quadrupole symmetry are Eq. H13(l and 
derivatives thereof; and on large scales we expect higher- 
derivative terms such as V^(V^ - V'^,2Vx^y)S\^ p] to 
be negligible. 

The primordial potential is related to the linear density 
field via: 



^p(k) = -47rG 



D{z) 



(14) 



where /b(w) = {wl. — w'^) / w^ . To get the -B-mode power 
spectrum, we replace Je with /_b(w) = 2'WxWy/w'^: 



2;^2 



^r(fc) 



bl /[/B(k2)-H/B(ki)] 






Note that the i?-mode intrinsic alignment power spec- 
trum contains no 0{Pg^^) contribution. This is because 
the tidal quadrupole field is a pure i?-mode (this is even 
true in the presence of nonlinear evolution), and produces 
only an £^-mode pattern of intrinsic alignments if a lin- 
ear model of galaxy ellipticities such as Eq. H13|l applies. 
Galaxy clustering modulates this field according to the 
galaxy distribution, and thereby transfers some power 
into _B-modes. 

The cross-power of the matter density and weighted 
shear is then found from Eq. (|15|l : 



Ps,^s{k) = -^Pt{k). 



B. Quadratic alignment model 



(18) 



The apparent ellipticity of a spiral galaxy is deter- 
mined principally by the orientation (and hence angu- 
lar momentum) of its disk. This angular momentum is 
believed to comes from external tidal fields perturbing 
the collapsing galaxy to form an anisotropic moment of 
inertia, which allows the galaxy to be "spun up" by a 
tidal quadrupole. In this case, the angular momentum 
vector L of the galaxy acquires an anisotropic probabil- 
ity distribution, and leads to a mean intrinsic ellipticity. 
This mean ellipticity vanishes to first order in the tidal 
field since a tidal field is required both to produce the 
anisotropic moment of inertia and then to apply a torque. 
The second-order contribution is :26] 



where p{z) is the mean density of the universe, D{z) oc 

(1 -I- z)D{z) is the re-scaled growth factor normalized to 

unity during matter domination, and G is the Newtonian J — C2 {T^^ — Ty 

gravitational constant. On linear scales, we have 5g = 

hgSun, S = Sur^, and hence the weighted intrinsic shear where the tidal tensor is 

is: 



IT T , 



.), 



7(k) = ^ 



Cj_p 
D 



{'^2x 



kly,2k2xk2y) 



K2 



5lin{^2 



+^<5.„(ki) 



d'ki. 



(15) 



T — 



47rG 



VpV^ - -(5^^V' 



S\^i 



(19) 



(20) 



where k2 = k — ki and we have chosen the wave vector 
k to lie on the a;-axis. (We are only interested in modes 
with k perpendicular to the line of sight.) The power 
spectrum of 7 is: 



(Since the tidal shears T are spin 2, we recall from the 
theory of addition of angular momenta that Eq.^|is the 
most general quadratic function of T with spin 2.) The 
resulting density-weighted intrinsic alignment is 



7Mk) = 



EE 



^1' 



(fc) 



x/i5(k2 



[Pt{k) 






'iki)Pi" 



(27r) 



[/£(k2) + /B(ki)] 

ki}, (16) 



C2P 

{2nfD 



hE{K,-k'2)SHn{K)SHn{K) 
b 



{2nY 



-SunO^z) 



d-'k'i d^k^(21) 



'(fc2),3 



where we have retained the convention that k is chosen to 
lie on the x axis, set kj, = 'k.a/\ka\ and kg = k — k'j^ — k2. 



and defined the /i-function by He — h^x — hyy with 



hx^{u, v) = ( u^u^ - -6^,^ j i vxVu - ^Sx„ ) . (22) 



The B-mode intrinsic aUgnment is obtained by replacing 
hE with Hb = 2hj:y (See Ref. 45] for a similar calcula- 
tion without the density- weighting 1 -I- i5g.) 

The predicted intrinsic alignment power spectra P^/ (fc) 
from this model can be derived by the same methods as 
used in Sec. IIII Al although since 7^ is second-order in 
^p (and hence Sun) there are more Feynman diagrams 
that contribute. The result is 



P^^^ik)^^{2 [hE{k,M)] 



+hlJ[hE{K, 



(27r)3 ^ 

hE{K,K) 



+ hE{KK)\' 

^ pj»(fci)pi™(fci)PJ-(fc^) . . . . 



(27r) 



d^k'id^k^j. (23) 



In Eq. H23|l we have defined k2 = k — ki and kg = k — 

k' — W 

K]^ K2. 

In the simple version of the quadratic-alignment model 
presented here, there is no density-shear correlation, i.e. 
5un and (1-1- 5g)^^ have zero cross-correlation assuming 
Gaussian 5un, linear biasing 5g = bgSun, and linear evolu- 
tion of the density field. We can see this as follows: since 
the three-point functions of Sun vanish, this correlation 
is equal to that of Sun and Sg-f^ . Of the contributing 
diagrams, those connecting Sun to the Sg vanish because 
(7^) — 0, hence we only consider the diagrams connecting 
6iin to one of the two tidal quadrupoles T^^ that make 
up 7^. But these vanish by symmetry since they must 
connect 6g, a scalar, to another tidal quadrupole T^i, eval- 
uated at the same point. This fact was first noted in the 
context of intrinsic alignment contamination of galaxy- 
galaxy lensing by Ref. 44j, which also pointed out that 
there can be nonzero cross-power Pgx^i{k) if we consider 
nonlinear density evolution. We will not consider these 
more complicated models here. 



C. Normalization 

We have normalized the intrinsic alignment amplitude 
Ci for the linear alignment model to match the shear 
autopower observed by Ref. [23 in the SuperCOSMOS 
data J43I. Ref. |23 gives the ellipticity variance a^{0) in 
square cells of side length 9 instead of the shear power 
spectra; these are related by 



'^m-n'Jji[^)fo[^ 



ic\ 



EE 



c, 



BB 



(24) 



where joiu) — sinu/u and 'R, = dejdrj « 2 is the shear- 
to-ellipticity conversion factor. We have assumed a pa- 
rameterized redshift distribution of the form 



dn 
dz 



oc z^e-fi-^^/^'") 



(25) 



where the median redshift for SuperCOSMOS is Zm = 0.1 
|29|. The SuperCOSMOS ellipticity variances are an or- 
der of magnitude greater than the weak lensing predic- 
tion and hence must be dominated by intrinsic align- 
ments and/or observational systematics; we have as- 
sumed here that the intrinsic alignments dominate, which 
provides an upper limit if observational systematics are 
not negligible. We have normalized to the largest-scale 
observation given by Ref. [23 , where Eq. l(T^ is best mo- 
tivated, and used the ellipticity covariance between the 
bj and R bandpasses as this is less sensitive to system- 
atics than the variances; this yields a^ = 1.6 x 10~^ at 
9 = 93.25 arcmin. 

For the linear intrinsic alignment model there is a dis- 
crete degeneracy in Ci since the sign cannot be deter- 
mined from the ellipticity variance alone; we have used 
in the plots the positive normalization Ci > because 
of its physical motivation (stretching of the galaxy along 
the tidal field). 



D. Results 

In Fig. 12 we have shown the estimated intrinsic align- 
ment contamination of the shear power spectrum for the 
linear alignment model. Panels (a-c) shows the contri- 
butions to the _B-mode shear power spectrum assuming 
a source redshift distribution given by Eq. H25|l with 
median redshifts Zm = 0.1, 0.5, and 1.0. In panels 
(d), (e), and (f) we have shown the intrinsic alignment 
contributions to the cross-spectra between different red- 
shift bins; only the gravitational-intrinsic correlation is 
nonzero since the redshift distributions do not overlap. 
For these panels we take the redshift distributions to be 
Gaussian with a^ — 0.1 (which is much less than the 
separation of the bins) . The contamination C^^ is only 
weakly dependent on CTz because the integral in Eq. Ullfl 
receives a contribution from a range of radii Ax oc az, 
whereas the integrand contains a factor of //3(x)j whose 
peak value scales oc crj^. 

The most striking result from Fig. [5] is that the 
gravitational-intrinsic contribution dominates the con- 
tamination of the shear power spectrum for the broad 
redshift distribution (Fig. O;). Even if the intrinsic- 
intrinsic shear correlations are small (in this case they 
are 0.5% of the lensing signal for the broad redshift dis- 
tribution with Zm = 1-0 at i — 500), the gravitational- 
intrinsic contamination to the cross-spectra can still be 
large (~ 5% at ^ = 500). This contamination is es- 
pecially pronounced for the cross-spectrum between the 
widely separated bins at z = 0.5 and 2.0 (Fig. Et): here 
|Cf VCf ^1 can be as large as - 30%. 



Shear power spectra for linear alignment model 
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(a) Power spectrum (Zf„=0.1) 



EE.total 

BB 

EE,GG 

EE,il 

EE,GI 



(b) Power spectrum (z^=0.5) 



(c) Power spectrum (z^=1 .0) 
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(d) Cross-power (z=0.5, 1 .0) 



(e) Cross-power (z=0.5, 2.0) 



(f) Cross-power (z=1 .0, 2.0) 
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FIG. 2: The gravitational {GG), intrinsic-alignment (//), and gravitational- intrinsic correlation {GI) contributions to the shear 
power spectrum, for the linear alignment model. The B-mode shear power spectrum is labeled "B" and has only an intrinsic- 
intrinsic contribution. The cross-term (GI) is negative for this model, so we have plotted its absolute values on the log scale. 
We show results for (a) the shear power spectrum for a survey with redshift distribution of Eq. H25|l with Zm = 0.1; (b) and 
(c) similar for for Zm ~ 0.5 and 1.0, respectively; (d) the shear cross-power between redshifts z — 0.5 and z = 1.0 (the slices 
have Gaussian redshift distributions with width cr^ = 0.1); (e) the shear cross-power between redshifts z = 0.5 and z — 2.0; 
and (f) the shear cross-power between redshifts z — 1.0 and z = 2.0. The intrinsic alignment amplitude Ci is normalized to 
SuperCOSMOS. Panel (a) is cut off at £ = 300, roughly the smoothing scale used for the intrinsic alignment calculation, since 
the model does not make sense on smaller scales. 



As noted above, there is no correlation between the 
gravitational and intrinsic shears in the quadratic align- 
ment model. Thus there is only intrinsic alignment con- 
tamination if the redshift distributions overlap, and the 
contamination is much lower because it scales as the 
square of the intrinsic alignment amplitude ex C2 instead 
of linearly as in the case of the GI cross-correlation. 



IV. METHODS TO ASSESS THE 
CONTAMINATION 



sic alignments, and shown that the interference between 
gravitational and intrinsic shears can dominate the con- 
tamination and reach levels that are important for high 
precision studies of weak lensing. In this section we dis- 
cuss two methods to estimate the amplitude of the ef- 
fect and separate it from weak lensing effect. The first 
method is based on density-shear correlations, while the 
second method is purely geometrical and uses only scal- 
ing of the cross-correlation signal with redshift. 



So far we have considered the contamination of the 
power spectra in cosmic shear surveys due to intrin- 



A. Density-shear correlations 

We are assuming there is a correlation between the 
intrinsic ahgnment and gravitational potential field that 
gives rise to shear: the latter can be reconstructed from 
the density field. On large scales we can assume galaxy 
density field is linearly proportional to the density field 
up to a constant, the so called bias parameter. In fact, 
there is no need to do the reconstruction itself: at the 
two-point function level the full information on this corre- 
lation is encoded in the density-shear correlation. We can 
thus constrain the intrinsic alignment models by compar- 
ing them to the observed density-shear correlation from 
galaxy surveys. 

To date these measurements have been done on small 
scales < lh~^ Mpc for the purpose of understanding con- 
tamination of galaxy-galaxy lensing. The quantity mea- 
sured in these studies is A7(r), defined as the mean tan- 
gential shear of a "satellite" galaxy a transverse distance 
r from the "primary." This is equal to 



Intrinsic shear correlation for linear model 



Aj{r) = 



/([I- 



s.mm 



^gix)hiix))dx\\ 



J{[l + 5g{o)]Sg{x))dx^^ 



(26) 



cf. Eq. (A5) of Ref. [ig. The idea behind this equa- 
tion is that we compute the galaxy-shear correlations for 
all close angular pairs on the sky and then use the cor- 
responding galaxy-galaxy correlation to estimate what 
fraction of these pairs is physically close also in three- 
dimensional space. The signal is then boosted by divid- 
ing by this fraction. 



Noting that (1 + Sg)jl 



7i, 



and taking the Fourier 



transform in the transverse dimensions, we get 



A7(r) 



I Ps^^^i{k)J2{kr)kdk _ JPs^^i{k)J2{kr)kdk 



JPs {k)Jo{kr)kdk bg J Ps{k)Jo{kr)kdk' 

(27) 
In Fig. Owe have plotted A7(r) as measured in the Sloan 
Digital Sky Survey (SDSS) [4^; the sample was domi- 
nated by objects at z « 0.1 and so we have displayed the 
model prediction at that redshift. (The data points were 
presented in Ref. 47] in physical separation units; we 
have converted them to comoving separation a,t z — 0.1.) 
We see from Fig. O that the model predicts a constant 
value of density-shear correlations on large scales, as ex- 
pected in a linear bias model. On scales smaller than the 
initial size of a galactic halo one does not expect much 
correlation. The transition scale is rather uncertain and 
depends on the somewhat arbitrary value of the cutoff 
we used in the calculation. A typical value should be 
given by the scale length within which the enclosed mass 
is of the order of a typical galactic halo mass, of order 
IO^^Mq. This gives fccutoff ~ (l-2)/i/Mpc, which are 
the two values used in Fig. El We see that the current 
constraints are inconclusive because they do not extend 
to sufficiently large scales. The results from this figure 
also suggest that galaxy-shear correlations are a powerful 
probe of identifying this contamination and that there is 
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FIG. 3: The intrinsic shear statistic A-f{r) of Eq. ^. The 
SDSS data points are shown, with the horizontal error bars in- 
dicating the range of radii used, the thick vertical error bars 
indicating la statistical errors, and the thin vertical error 
bars representing 99.9% confidence limits including systemat- 
ics (principally shear calibration and removal of lensing signal) 

na. 



plenty of statistical power in the SDSS data to estimate 
the contamination by extending the analysis of Ref. |43 
to larger transverse separations. 



B. Geometric projection 

The analysis of intrinsic alignment contamination 
should not be based solely on low-redshift information: 
the local universe may not accurately reflect the intrin- 
sic alignments of the weak lensing source galaxies since 
(for example) mergers could cause the intrinsic align- 
ment signal to evolve with redshift. Consistency of shear 
cross-spectra obtained from different types of galaxies 
in the same range of redshifts would argue against sig- 
nificant gravitational- intrinsic correlations, but if an in- 
consistency is observed that is not attributable to ob- 
servational systematics, then some method of separating 
pure lensing {GG) from gravitational-intrinsic {GI) sig- 
nals will be needed. 

One possibility is to take advantage of the different red- 
shift dependence in the GG and GI signals when applied 
to tomography between slices at redshifts Za < zp (the 
redshift slices should be separated far enough to elimi- 
nate overlaps in their distribution and hence // contam- 
inaiton) . In the idealized case where the photo-zs for the 
nearby redshift bin Za have negligible uncertainty, the 
GI signal at a fixed i and Zq, rises with increasing zjs as 






/?)cx 



Da{0,z) 



[i + C(^)] 



dn 



I dz, 
dzj p 



(28) 

where Da is the comoving angular diameter distance sep- 
arating two redshifts and C(z) is the shear calibration er- 
ror for galaxies at redshift z (which ideally would be made 
negligible). The GG (pure lensing) signal also rises with 
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increasing zp, but more slowly because z^ in Eq. 1^5)1 is 
replaced by the redshift of the contributing structures, 
which must be < Zq,. In principle, one could project out 
signals in tomography that have the redshift dependence 
of Eq. (EHll- 

If the effective ratio of angular diameter distances gafi 
is measured by cosmography - i.e. measurement of the 
relative shear signal at different zp behind a sample of 
lenses at Za - then provided that the same source sam- 
ple at Zj3 is used for tomography as for cosmography, the 
right-hand-side of Eq. if^ is obtained directly ^3j Ull 
up to an irrelevant multiplicative factor. This geometri- 
cal method has the practical advantage of automatically 
accounting for any shear calibration errors, even if they 
are redshift-dependent, and eliminating the need for ex- 
ternal information about the angular diameter distance 
function Da- Since the density-intrinsic shear correla- 
tion may evolve with redshift, the de-projection of signals 
proportional to Eq. if^ will have to be done separately 
for each redshift slice z^- The method is more accurate 
if the slices are narrow so that the effective Za for the 
intrinsic alignment does not differ significantly from the 
effective Za for cosmography, thus precise photo-zs would 
be extremely valuable. 

The geometric method is model-independent in the 
sense that it works for arbitrary k and z-dependence 
of the density-intrinsic shear correlation, but of course 
some cosmological information is lost. An alternative ap- 
proach would be to marginalize over a parameterization 
of Ps^^! (fc, z) (perhaps forced to match on to the Ps.^' {k) 
measured in the local universe at z w 0); this may retain 
more cosmological information if Ps^^i {k, z) is e.g. forced 
to be a smooth function of z, but of course it sacrifices 
the model independence of the geometric technique. 



V. DISCUSSION 

In this paper, we have considered the contamination 
of the power spectrum in cosmic shear surveys due to in- 
trinsic alignments, and shown that the gravitational and 
intrinsic shears need not be independent. For some mod- 
els, the contamination is dominated by the gravitational- 
intrinsic cross-power rather than thet power spectrum of 
the intrinsic alignments themselves. One implication of 
this is that the B-mode signal, frequently used as a sys- 
tematics test, may be a misleading indicator of intrin- 
sic alignment contamination, since it is sensitive only to 
the intrinsic-intrinsic autopower C^ ' but not to the 
gravitational-intrinsic cross-power Cf^^ , which can only 
be nonvanishing for the E modes. Fig. |2Ic) provides an 
extreme example: here the change in Cf^ due to in- 
trinsic alignments exceeds the i?-mode power spectrum 
Cf^ by > 2 orders of magnitude. The opposite extreme 
is the quadratic alignment model in which there is no 
gravitational-intrinsic correlation and the E and B mode 
contributions from intrinsic alignments are similar. 

A realistic model of intrinsic alignments will likely fall 



somewhere in between these two extremes. The lin- 
ear alignment model (Eq. I13|l obviously predicts maxi- 
mal correlation between the density and intrinsic shear 
fields on large scales in the sense that the correlation 
coefficient ps^Xfi{k) = Ps;yi{k)/y/Ps{k)P;yi{k) has abso- 
lute value w 1. Any higher-order corrections to Eq. H13|) 
will therefore reduce the correlation coefficient. On the 
other hand, the simplest form of the quadratic alignment 
model (Eq. I19|) is probably too optimistic, even though 
the model is well-motivated by tidal torque theory. We 
found Ps;yi{k) = for this model not because of any fun- 
damental symmetry principle, but rather because the as- 
sumed Gaussianity of the density field causes third-order 
statistics in 6 to vanish. In the real universe we can ex- 
pect tidally torqued galaxies to have Pg^ji (fc) -# because 
of non-Gaussianity from nonlinear evolution |4J| , and be- 
cause of third- and higher-order corrections to Eq. H19() 
that we have neglected. 

While we have presented these results for the power 
spectrum, the non-Gaussianity of the density field, and 
hence of cosmic shear, has motivated studies of other 
statistics, most notably the bispectrum J53, ^M, HE 
153 . I54L l55| . The bispectrum can in principle pick up 
three types of intrinsic alignment contamination: the in- 
trinsic shear bispectrum {III), and the gravitational- 
gravitation-intrinsic (GGI) and gravitational-intrinsic- 
intrinsic (Gil) cross-bispectra. If tomography is used, 
the measured bispectrum components are Bf^f'g (a/37), 
where a, (3, and 7 are indices indicating the source bins 
centered at redshifts Zq < zp < Zj. Assuming that 
these bins do not overlap significantly in redshift, the 
/// contaminant is only nonvanishing for the triplets 
where all three source galaxies lie in the same redshift 
slice Za — zp — Zy] this is only a small fraction of the 
triplets and hence little statistical power is lost by reject- 
ing them .33]. GGI is produced by a density-density- 
intrinsic shear bispectrum in the nearest of the three 
slices, and it can in principle contaminate any tomo- 
graphic shear bispectrum in which Za < zp < z^y. Gil 
is produced by the density-intrinsic-intrinsic bispectrum, 
again in the nearest of the three slices, but this time in 
order to have the intrinsic shear be correlated and to pro- 
duce a lensing effect on the more distant slice, Gil sat- 
isfies a "selection rule" Za = zp < z^. Thus every triplet 
of source screens is potentially contaminated by exactly 
one of GGI, Gil, or III. Note also that the quadratic 

which at lowest order pre- 



alignmcnt model (Sec. IIII B|l 
dieted no GI contamination to the power spectrum, does 
predict a GGI contamination Bggyyi{ki,k2,k3) ^ be- 
cause the intrinsic shear 7^ is a quadratic function of 6. 
We defer a calculation of the magnitude of this effect to 
a future paper. 

Our analysis clearly motivates observations to con- 
strain the density-intrinsic shear correlation. Just as 
the intrinsic shear autopower can be constrained by low- 
redshift measurements where the lensing signal is negli- 
gible 24, 29], so it should be possible to use these same 
low-redshift measurements to measure Pg^^yi (fc) provided 



that a tracer of the density field is available. The num- 
ber density of galaxies is one possibility; recent observa- 
tions have probed the Kalaxy-intrinsic shear correlation 
on 0.02-lh-^ Mpc scales [mIs^, and it would be useful 
to extend these studies to larger scales where linear bias- 
ing is valid. (At least one density-shear correlation mea- 
surement is available on large scales [23 • These authors 
were interested primarily in using intrinsic alignments 
to reconstruct the density field; they therefore did not 
use lensing shear estimators, and they measured three- 
dimensional separations in redshift space, so their results 
are difficult to interpret in the present context.) For such 
observations, it would also be useful to distinguish types 
of galaxies (e.g. early vs. late); if one type of galaxy turns 
out to have significantly greater Ps^^i (k), the shear power 
spectra could be performed principally on the other types 
so as to reduce contamination. 

A second approach to eliminate the contamination 
is to project out the information that correlates with 
the known redshift scaling of the signal in the cross- 
correlations between narrow redshift bins, assuming these 
are available from the photo-z methods. This method is 
geometric and achieves perfect decontamination under 
these idealized assumptions, but projection does destroy 
some information in the data. It also puts severe require- 
ments on the accuracy of photometric redshifts, which 
may require many passbands to achieve this 57]. 

In summary, we have shown that cosmic shear surveys 
will need to consider the possibility of intrinsic alignment 
contamination from its interference with weak lensing in- 
duced shear. The magnitude of the effect is extremely 



model-dependent and can range from essentially zero 
(simplest quadratic alignment models) to severe (sim- 
plest linear alignment models). In the latter case it can 
exceed intrinsic correlations by an order of magnitude 
and cannot be identified using _B-mode power spectrum, 
nor can it be eliminated using cross-correlations between 
different redshift slices. This result is disappointing and 
implies that this effect could be very damaging for efforts 
to use weak lensing as a high precision test of cosmol- 
ogy with future surveys such as Pan-STARRS, LSST or 
SNAP. However, the situation is far from hopeless: here 
we have considered several possible ways to constrain 
and/or suppress the contamination. The next step obser- 
vationally is a density-shear correlation analysis with the 
shear computed using a lensing estimator in a wide-angle, 
shallow survey such as the Sloan Digital Sky Survey 53] 
to constrain Pg;yi{k,z <^ 1); the next step theoretically 
is a thorough analysis of the various methods described 
above for suppressing the gravitational-intrinsic contam- 
ination, in particular to understand how much degrada- 
tion on the cosmological parameters occurs and whether 
they increase sensitivity to any of the other potential sys- 
tematics such as photo- z errors, shear calibration errors, 
and spurious power. 
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